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Interplay between single-particle and two-particle tunneling in normal 
metal-fi-wave superconductor junctions probed by shot noise 
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We discuss how life-time broadening of quasiparticle states influences single- and two-particle current transport 
through zero-energy states at normal metal/d-wave superconductor junctions. We distinguish between intrinsic 
broadening (imaginary part 77 of the energy), which couples the bound states with the superconducting reservoir, 
and broadening due to leakage through the junction barrier, which couples the bound states with the normal metal 
reservoir. We show that shot noise is highly sensitive to the mechanism of broadening, while the conductance is 
not. In the limit of small but finite intrinsic broadening, compared to the junction transparency D, tj/Aq <C D, 
the low-voltage shot noise at zero frequency and zero temperature becomes proportional to the magnitude 77 of 
intrinsic broadening (Ao is the maximum d-wave gap). 
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1. Introduction 

Tunneling experiments have for a long time 
been an important probe of superconducting 
properties According to early theoretical 

analyses in terms of the tunneling Hamiltonian 
model[^, the conductance of a normal metal- s- 
wave superconductor tunnel junction is propor- 
tional to the superconductor density of states. 
Within this model only single-particle tunneling 
is taken into account. For more transparent inter- 
faces. Blonder, Tinkham and KlapwijkQ] (BTK) 
presented a way to calculate the conductance 
within Bogoliubov-de Gennes quantum mechan- 
ics. Their approach is valid for arbitrary trans- 
parency D of the interface, but assumes that the 
junction system has no impurities and that inelas- 
tic scattering is negligible. BTK showed that for 
intermediate-to-high transparency, so-called An- 
dreev reflection Q| plays an important role in de- 
termining the conductance. Andreev reflection 
is the process where an electron with an energy 
near the Fermi surface, incident on a supercon- 
ductor from a normal metal, is converted into a 
hole with essentially unchanged momentum, but 
with opposite group velocity. The reversed pro- 



cess is also possible. In the Andreev reflection 
process a Cooper pair is formed on the supercon- 
ductor side and a charge 2e is transferred over the 
NS interface, which corresponds to a two-particle 
tunneling event In the low-transparency limit, 
Andreev reflection is suppressed since it is of sec- 
ond order in the transparency D, while ordinary 
single particle tunneling is of first order in D. 
Therefore, in the tunnel limit, D < 1, the BTK 
result reduces to the tunnel model result. 

In contrast to the clear situation for s-wave 
junctions described above, the interpretation of 
the conductance of c?-wave NIS junctions|^-|j 
where interface resonant states are present is 
more complicated. Surface/interface states with 
zero energy are necessarily formed as a conse- 
quence of the d-wave symmetry of the order pa- 
rameter, if the interface is oriented so that the 
order parameter changes its sign along quasipar- 
ticle trajectories involving normal reflection at 
the interface!^. For a recent review of the im- 
plications of these surface/interface zero-energy 
states (ZES), see Ref. It is important that 

two physically distinct situations can be realized 
depending on the relation between (i) intrinsic 
broadening due to e.g. surface roughness, im- 
purities, or phonons, which connects the inter- 
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face states with the bulk superconductor, and (ii) 
broadening due to leakage to the normal reservoir 
through the barrier. When intrinsic broadening 
is negligible, i.e. when the interface states are de- 
coupled from the superconducting reservoir, the 
interface states can participate in current trans- 
port only via the Andreev reflection process. In 
fact, the Andreev current is resonantly enhanced 
by the ZES, and is of first order in the trans- 
parency D. On the other hand, when the intrinsic 
broadening is dominating, the interface states will 
instead assist single particle tunneling. Thus, for 
a d-wave superconductor with ZES, a zero-bias 
conductance peak (ZBCP) is present in both lim- 
its, but the type of tunneling responsible for it is 
very different. Even if both types of life times are 
large, it is the quotient 
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which will determine what type of tunneling is 
dominating: single particle tunneling if g ^ 1 and 
two-particle tunneling if q -C 1 [|lO|,0 ■ The intrin- 
sic relaxation time Tr (x h/rj is set by the damping 
introduced by an imaginary part r] of the quasi- 
particle energy, while the life time oc Ti/DAq 
due to leakage through the barrier is determined 
by the transparency D of the barrier (Aq is the 
maximum d-wave gap). 

In the present paper we will discuss the in- 
terplay between single-particle and two-particle 
tunneling through the ZES, with emphasis on 
the limit of small but finite damping. We will 
show that it is possible to discriminate between 
the two types of tunneling in a measurement of 
low-voltage shot noise at zero temperature and 
zero frequency. Shot noise was recently calcu- 
lated for c?-wave NIS junctions in Refs. [p^ , p^ , 
under the implicit assumption of negligible intrin- 
sic broadening {q = 0). It was found that the 
zero- voltage differential shot noise level is zero, 
{dS/dV){V = 0) = 0, because the effective An- 
dreev reflection probability is enhanced to unity 
near zero energy by the ZES resonance. Below 
we will study the effect of intrinsic broadening, 
and show that the noiseless character of the ZES 
is quickly lost when rj increases from zero. This 
happens because the effective Andreev reflection 



probability is reduced from unity with increas- 
ing broadening; single-particle tunneling via the 
ZES then becomes possible on the expense of two- 
particle tunneling, and fluctuations between the 
two channels is introduced. When broadening is 
small, the low-voltage shot noise will be propor- 
tional to rj, which makes it possible to probe the 
magnitude of intrinsic broadening by a measure- 
ment of noise. 

2. Calculation of current and shot noise 

In our calculation we use the coherent BTK 
scattering approach for a specular d-wave NIS 
junction, modified to include intrinsic broaden- 
ing on a phenomenological level by an imaginary 
part of the energy. We assume that the normal 
metal is at a:; < 0, while the superconductor is at 
a; > 0. We let a c?-wave gap node point towards 
the interface, so that the spectral weight of the 
ZES is large. It is known that for this orienta- 
tion the gap is substantially suppressed near the 
interface ||l4]. However, in order to make clear 
our points on the presence of an interplay be- 
tween single-particle and two-particle tunneling 
as a function of q, and that shot noise will be sen- 
sitive to the nature of tunneling, it is sufficient to 
assume a step function form of the gap function 
A{6,x) — 9(x)Aosin20, where 9 is the angle of 
quasiparticle propagation relative to the interface 
normal. The exact form of the specular barrier is 
not important and we assume a (5-function poten- 
tial characterized by a transmission probability 
D{6) = cos^ 6 / {cos^ 9 + Z"^), where the dimen- 
sionless quantity Z is a measure of the strength 
of the potential. 

In the original BTK theory, effective normal 
reflection Rn and Andreev reflection Ra proba- 
bilities are calculated and then the charge current 
spectral density is expressed through them. We 
will follow a slightly different route and ex- 
press the charge current density in terms of prob- 
ability current densities flowing along the energy 
axis. This viewpoint makes it possible to rigor- 
ously divide the total charge current I{V) into 
single-particle Ii{V) and two-particle hiV) cur- 
rents. When an electron at energy E is injected 
from the normal reservoir into the junction, it 
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Figure 1. Quasiclassical paths illustrating the 
structure of the scattering state in real space 
(the xy-plane) due to an incoming electron from 
the normal metal side (N). In between the bar- 
rier (I) and the superconductor (S), we introduce 
a small normal region (N). Normal scattering 
takes place at the insulator and Andreev reflec- 
tion takes place at the NS interface. 



is accelerated by the voltage drop eV over the 
barrier. At the superconductor it is Andreev re- 
flected into a hole which, having the opposite 
charge, is also accelerated by eV when it tun- 
nels over the barrier back to the normal metal. 
Thus, the Andreev reflected hole is at the en- 
ergy E + 2eV. Each particle carries a certain 
probability current, proportional to the absolute 
value squared of the wavefunction. Thus, it is 
natural to change focus from the charge current 
flow through the junction (along the a;-axis), to 
the probability current flow along the energy axis 
{E-axis), see Figs. ^^j. In the figures, we have 
introduced a small auxiliary normal region (N) 
in between the barrier (the insulator, I) and the 
superconductor (S). The final results will not be 
affected by the presence of N and in the end we let 
the thickness of this region go to zero. On each 
of the electron and hole legs in N (see Fig. ||), 
we can define the total probability current densi- 
ties and j^, respectively, including the effects 
of normal backscattering at the insulator. For 
energies outside the gap, the Andreev reflection 
probability is less than unity, which leads to a 



Figure 2. Structure of the scattering state in en- 
ergy space (the a;i?-plane). 



leakage — jf^ of probability current during An- 
dreev reflection. The total leakage is the single- 
particle current. The rest of the probability cur- 
rent (the part surviving Andreev reflection) con- 
tinues out into the normal metal as a hole current, 
and contributes to the two-particle current. For 
subgap energies, the Andreev reflection probabil- 
ity is unity; probability current is then conserved 
during Andreev reflection, the leakage jl—j^ van- 
ishes, and single particle tunneling is quenched. 

In the zero-temperature limit, the current then 
takes the form 115(1 



eRnI{V) 
h{V,9) 



eRn[h{V) + h{V)] 

1 /•'^/2 



2{D) 



■n/2 



decos6[h{v,e) + i2{v,6)] , 



-eV 

2 [ dE- 

J-eV 



|1 -i?aa|2 



where Rn — ■nh/ e^kpLy {D) is the normal state 
junction resistance, Ly is the junction width, kp 
is the Fermi wave vector, (D) = J dO cos9D{9)/2, 
and a — a{9, E) is the Andreev reflection ampli- 
tude. The amplitude d is calculated at the angle 
vr— 0, the angle of propagation after normal reflec- 
tion at the barrier. The differential conductance 
G{V) is obtained by differentiation with respect 
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Figure 3. (a) The total current, (b) the differ- 
ential conductance, (c) the single particle cur- 
rent, and (d) the two-particle current, for a d- 
wave NIS junction oriented so that a gap node 
is pointing towards the interface. The solid 
lines are calculated for 77 = 0, while the dashed 
lines are calculated for 77 — O.OOIAq. The 
zero-energy states resonantly enhance the two- 
particle current, which is of first order in the 
transparency. Increasing broadening enhances 
the single-particle current on expense of the two- 
particle current, but leaves the total current un- 
affected. The transparency of the junction is 
{D) = 0.026, and the temperature is zero. 



to V. 

For the calculation of shot noise we follow the 
technique in Ref. see also ||l2|. The formula 
for the noise, here limited to zero temperature 
and zero frequency, is 



RnS[V) 



1 



t/2 /-O 

, d0cos9 / dES{0,E), 

\D) J-Tr/2 J-eV 

RNil-RN) + RA{l-RA) (3) 

+ 2R]yRA, 



where Rn ~ R\i — ^ Raa\^, and Ra = 

iP^lap/ll — Rad\^. In the subgap region, for zero 
temperature and negligible broadening, the only 
surviving source of noise is a term involving scat- 
tering states separated in energy by 2eV: the first 



scattering state consists of an electron incoming 
from the normal metal reservoir at the energy E, 
Andreev reflected as a hole which is ejected back 
into the normal metal at the energy E + 2eV; the 
second scattering state consists of a hole incom- 
ing at E + 2eV, Andreev reflected as an electron 
emerging at E. 

We plot the total current I{V), the conduc- 
tance G{V), the single-particle current Ii{V), and 
two-particle current hiV) in Fig. ||(a)-(d). From 
Fig. ^ it is clear that in the absence of broad- 
ening, the solid lines {q = 0), the ZBCP is due 
to the two-particle current only. One can show 
that the two-particle current is of first order in 
the transparency, I2 oc D, because of the ZES 
resonance. The width of the ZES resonance is 
r(e') = D{e)\A{0)\/ h^/RiWj for each quasipar- 

ticle trajectory angle 9, see also e.g. Ref. [^Ojjl^. 

To study the effects of intrinsic broadening, the 
quasiclassical Green's function technique [p^ can 
be used. Recently, a highly useful parameter- 
ization was introduced by Nagato, Nagai, and 
Hara|19| for the equilibrium case, and by Es- 
chrig|20| for the non-equilibrium case, which can 
be used to separate electron and hole parts of the 
Green's functions. Probability currents flowing 
along the energy axis are then easily found, and 
one can divide the total charge current into single- 
and two-particle currents, in a one-to-one map- 
ping to the approach outlined above. Eventually, 
the result is the same expression as in Eq. (H), 
but with the Andreev reflection amplitudes sub- 
stituted by generalized Andreev reflection ampli- 
tudes containing the effects of damping in the 
form of an imaginary part of the energy [ pO| , pl| . 
When broadening increases, the two-particle cur- 
rent is suppressed and the single-particle cur- 
rent is enhanced. However, the total current and 
the conductance are not particularly affected, as 
shown in Fig. |^(a)-(d) by the dashed lines. In 
fact, in the limit of large broadening, q ^ 1, 
the ZBCP is solely due to single-particle tunnel- 
ing. In this case, it is convenient to instead of 
the above approach apply the tunnel formula for 
the conductance and calculate the local density 
of states at the surface, as in e.g. Ref. |22|. In 
this limit the width of the ZBCP is directly pro- 
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Figure 4. Differential shot noise dS/dV at zero 
temperature and zero frequency, calculated for 
the same junction as in Fig. ^. In (a) for ry = 0, 
and in (b) for r] = O.OOIAq. The noiseless charac- 
ter of the zero-energy states is quickly lost when 
broadening is introduced. 



portional to the magnitude of the broadening rj. 

In Fig. ^ we plot the differential shot noise, 
dS/dV . The ZES resonance does not produce 
any noise [dS/dV){V = 0) = 0, see Fig. |(a), 
as found in Refs. [ p^[p^ ]. However, when a 
small amount of intrinsic broadening is intro- 
duced, the noiseless character of the ZES is lost, 
see Fig. ^b). This behavior can be understood 
in the following qualitative way. When there is 
no broadening, the ZES resonantly enhances the 
Andreev current. The effective Andreev reflection 
probability is unity, despite the fact that D ^ 1, 
c.f. Eq. d) ioYE = and Fig. |(d). The result of 
zero noise follows directly from Eq. (||) , since the 
relation Rn -\- Ra = 1 holds in the subgap region 
(probability current conservation). When broad- 
ening increases from zero, the effective Andreev 
reflection probability is reduced from unity, and 
the noise becomes finite. In the limit of small 
intrinsic broadening, q <C 1, the noise at zero- 
voltage is proportional to the magnitude of broad- 
ening, {dS/dV){e,V = 0) oc r7/r(e'), where T{0) 
is the width of the ZES in the absence of damping. 



3. Summary 

We have discussed the interplay between single- 
particle and two-particle tunneling through zero- 
energy states (ZES) in d-wave NIS junctions. For 
small intrinsic broadening of the ZES, compared 
to the broadening due to leakage over the barrier 
to the normal metal reservoir, q <C 1, current is 
transported via two-particle tunneling. On the 
other hand, for large damping, q ^ 1, only single 
particle tunneling is present. We have shown that 
shot noise is highly sensitive to the type of tunnel- 
ing, although the conductance is not. The noise- 
less character of the ZES found in Ref. |p|,|3| 
is quickly lost when damping is introduced. For 
small but finite intrinsic broadening, q <IC 1, the 
low- voltage shot noise is directly proportional to 
the magnitude of broadening. 
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